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Abstract 

The introduction of a modulus z{K), analogous to u = {trcj)"^) in the N = 2 SUSY SU(2) 
gauge theory solved by Seiberg and Witten, and whose defining property is the invariance 
under the symmetry and duality transformations of the effective coupling K, reveals an 
intriguing correspondence between the D = 2 Ising and Potts models on the square lattice. 
The moduli spaces of both models, the spaces of inequivalent effective temperatures K, 
correspond to a three-punctured sphere Ai^ = P^(C)\{z = ±l,oo}. Furthermore, in both 
models, the locus of Fisher zeroes is given by the segment joining Zc = —1 to Zc = +1. 
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1 Introduction 



Duality plays a crucial role in the latest developments of string theory and has proven to 
be a tremendously powerful tool in the solution of non trivial D = 4 supersymmetric field 
theories, namely the Seiberg-Witten models [|I|. 

In 0, Fendley and Saleur showed that various (1 + 1) -dimensional quantum impurity 
systems display an exact form of self-duality, analogous to the famous Kramers- Wannier 
duality relation in the D = 2 Ising model They achieved this by expressing the relevant 
quantities in terms of contour integrals over certain hyperelliptic curves. 

Previously, in 0], Fendley had expressed the magnetization in the Kondo model and the 
current in the boundary sine-Gordon model in terms of such contour integrals, noticing the 
similarity with the Seiberg-Witten results. 

Fendley and Saleur also posed the question whether the duality relation itself is enough 
to solve the specific model, without making use of the Bethe ansatz technique. 

This problem motivates our analysis which will focus on the D = 2 Ising and Potts 
models on a square lattice. We will see that the Kramers-Wannier duality relation and 
its generalization to the Potts model naturally lead to identify a three-punctured Riemann 
sphere Ai^ = F^{C)\{z = ±1, oo} as the moduh space of the given model. 

This is achieved by the introduction of a modulus or uniformizing coordinate z, with the 
property of being invariant under the symmetry and duality transformations. 
In particular 

_ 1 + sinh^(2fs:) 
^ ~ 2sinh2(2ir) ' 

z{K) = z{K*) = z{-K) = z{K + z^), (1) 

in the Ising case, and 

' - 2^q{e^ - 1) ' 

z{K) = z{K*) = z{K + 2m), (2) 

in the Potts case, where K is the effective coupling between neighbouring spins. 

Furthermore, in terms of the appropriate modulus z, the locus of Fisher zeroes both 
in the Ising and in the g-state Potts model is given by the segment joining the complex 
temperature singularity z = —1 and the physical singularity z = +1. 
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2 Ising model 

The Ising model is defined as follows [Q. At each site of a /^-dimensional lattice, there is a 
spin variable cxj = ±1 interacting with its nearest neighbours only through the energy 

H = -J^cTiaj, (3) 

where the symbol denotes the sum over nearest neighbours pairs. The partition function 
Z is given by 

z=t.^-'''=t.-w(kj:^.^], (4) 

W} W} \ {id) J 

where K = j3J and the sum is over all possible spin configurations. 
The free energy per spin f{K) is given by 

-/3/W=mn!^, (5) 

Let us define 

F{K) ^ -Pf{K). (6) 

In the following, we will consider a two-dimensional square lattice with sites. 
At high temperature and correspondingly small K 

Zugh{K) = E n cosh(i^)(l + a^a^ tanh(ir)). (7) 

The sum over all ai = ±1 selects only those terms with even powers of cTj, while the others 
cancel exactly. We can represent each term diagrammatically by a line connecting the sites 
i and j for each factor tanh(A')(Ti(Tj. Therefore, the non vanishing contributions come from 
all closed loops on the lattice 

Zf,,g,{K) = 2^(cosh Yl (tanh Ky. (8) 

loops 

where / denotes the length of the loop, that is the number of bonds in the loop. Note that 
the loops can be disconnected. By (H), we see that the partition function is invariant under 

K ^ -K, (9) 

due to the fact that / is even. 
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Likewise, at low temperature, large K, a spin configuration is identified by the boundaries 
of positive spin droplets in a negative spin background or viceversa. These boundaries form 
loops and the partition function can be expressed as 

ZUK) = 2e^^^ E e-^^^ (10) 

loops 

From this we easily see that under the shift 

+ (11) 

the free energy transforms like 

Fi,^{K + i-) = FiUK) + Y- (12) 
Furthermore, by (||) and ([lO|) , we see that if we set 

= tanh (13) 

the high-temperature and low-temperature expansions are mapped into each other 

(14) 



2e(2iVi^*) - 2^cosh(ir)2JV 

This is the Kramers- Wannier self-duality relation 0. In terms of the free energy F{K), 
(|14D becomes 

Fio^iK*) - ^lnsinh(2K*) = F,,igh{K) - ^ lnsinh(2if), (15) 

where we used 

= tanh K < — ^ sinh(2i^*) sinh(2ir) = 1. (16) 

Arguing that F{K) could have only one singularity for > 0, Kramers and Wannier 
concluded that the transition temperature had to coincide with the fixed point K*{Kc) = 
under the duality mapping 

sinh2(2irj = 1 ^ = ^ In (1 + ^2). (17) 



Note that the duality mapping (|T6[) connects the low temperature region with the high 
temperature region and viceversa: in this sense it is a mapping between a strong coupling 
regime and a weak coupling one. Furthermore, it is involutive 

K*{K*{K)) = K. (18) 
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3 The modulus z 



In Seiberg and Witten managed to calculate the low-energy Wilsonian effective action of 
the N = 2 supersymmetric SU (2) Yang-Mills theory. This effective action is given in terms 
of the so-called prepotential J-'{a), which is a function of the vacuum expectation value a 
of the scalar field contained in the supersymmetric multiplet. Most importantly, N = 2 
supersymmetry constrains jF(a) to be a polymorphic function, that is a multivalued analytic 
function. 

The crucial feature is that the theory has a moduli space Ai, namely a manifold of 
physically inequivalent vacua, which is parametrized by the gauge-invariant coordinate 

M=(tr02). (19) 

The low-energy description breaks down at certain points Uc, due to the appearance of extra 
massless particles. Furthermore, as u loops around these singular points, a{u) and its dual 
o-niu) = undergo a monodromy, they transform into a linear combination of themselves 



a(u) 




(20) 



where u = Uc + e'^^^{u — Uc) and A, B, C, D are integers satisfying AD — BC = 1. 

These monodromies correspond to duality transformations relating two theories with 
different effective coupling constants r = 9^jF(a), Imr > 0.. Eq. (|20| ) implies that 



r 



daD At + B 



and 

u{^t) = u{t). (22) 

Arguing that the theory had only three singularities, which could be set at {uc = ±1, oo}, 
Seiberg and Witten were able to determine a{u) and a£,{u) exactly in terms of period integrals 
of a suitable meromorphic 1-form over the two canonical homology basis cycles of an elliptic 
curve parametrized by u. 

These period integrals have the desired monodromies and are known to satisfy second- 
order linear differential equations called Picard-Fuchs equations. In particular 




d'u + T7-^ ?/ =0. (23) 



Note that, when u loops around one of the punctures, the differential equation does not 
change because the potential is single-valued. Therefore, a^iu) and a{u) must be a linear 
combination of aoiu) and a{u) with constant coefficients, as in (pOl). In general, multival- 
ued functions with non-trivial monodromies are naturally associated to linear differential 
equations with meromorphic potentials. 

From Eq. (|23|) , Matone derived a non-perturbative identity relating u and JF(a) 

u = Txi (j^{a) aao] , (24) 
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which allowed to find recursion relations for the instanton contributions to the prepotential 
The identity (24) was verified by instanton calculations [12] [13], and it follows from the 



superconformal Ward identities as well |TJ]. In Bonelli, Matone and Tonin took it as a 
starting point for a rigorous derivation of the Seiberg-Witten result by reflection symmetry, 
without any assumption on the number of singularities |T3[. They univocally identified the 
three-punctured sphere Ais = P^(C)\{mc = ±1, oo} as the moduli space of the model. 

The Kramers-Wannier self-duality relation and its generalization in the Potts model 
imply that the physical properties of the system at low and high temperatures, equivalently 
high and low K, are related. 

Furthermore, K plays the same role as the effective coupling constant r in the Seiberg- 
Witten model. Therefore, it is natural to look for a new variable z{K), analogous to u{t), 
that is left invariant under both the duality mapping (p!6D and the symmetries (|^), (0) 

z{K)=z{K*)=z{K+—) = z{-K). (25) 

Consequently, we may regard the free energy F{K) as a multivalued function of z, J-'{z), and 
the various transformation properties of F{K) under (P), (|ll|) and (|16|) will be reflected in 
the polymorphicity of T{z). In fact, performing one of these mappings is equivalent to the 
modulus z going around a non trivial closed path in the moduli space M.. Correspondingly, 
the free energy J-'{z) will undergo a monodromy. 

This procedure will show us that the natural setting or moduli space for the Ising model 
in two dimensions is again a three-punctured sphere M.^ = F^{C)\{z = ±1, oo}. As we will 
discuss below, this is of course encoded in Onsager's solution 0]. We will consider the Potts 
case later. 

By (^), it is natural to consider the effect of the various transformations on s = sinh(2i^'). 
In particular, s —s under both (|) and (|ll]), and s — > under (|l^). We see that the 
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maps 

f : s ^ -s g : (26) 

are involutive and commute 

fof=id = gog, f o g = g o f ^ {f o g) o {f o g) = id. (27) 
Hence the polynomial in the variable x 

Fix, s) = ix- s)ix - fis))ix - g{s)){x - /(^(s))), (28) 

will satisfy 

P{x, s) = P{x, f{s)) = P{x, g{s)) = P{x, f{g{s))), (29) 
and its coefficients will be invariant under (||), (0) and (|16|). We have 

P{x,s) = {x-s){x + s){x--){x + -) = x'- ( +1, (30) 

s s V / 



which implies that is the quantity we are looking for. We remark that there are no 
transformations other than (^), ( pJ]) and (|16D which leave z unchanged. This is a minimal 



choice in a sense, since no extra spurious transformations enter the picture. For instance we 
can easily check that all the quantities: 

His,n) = ^, (31) 

are invariant under (P), (|ll|) and (^). But these are not the only symmetry transformations. 
For example, H{s, 2) would be invariant under s — > is as well: however, we do not know 
how the free energy transforms under this map. 
Finally, we will set: 



- = (32) 



+ 1 

so that the critical temperature corresponds to Zc{K) = 1. 

We remark that there is no loss of generality in picking s instead of tanh(ir) or e^^^ as 
the building block. The invariants one finds using the latter are equivalent to z. 

We would like to view the free energy as a function of z. To this end, we have to invert 
( p2D and obtain s = s{z), which is of course a multivalued function of z. In general, by 
construction, for a given value of z there are four different values of s, obtained by solving 
the equation 

A{s, z) = s'^- 2zs'^ + 1 = 0. (33) 
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However, there are certain critical values Zc such that the polynomial (Q) has multiple roots. 
Equivalently, z^. = z{si), where Sj are ramification points of the map z{s). As z winds around 
one of these critical values, we move from one branch of s{z) to another. This is analogous 
to = x: as X winds around or equivalently oo, we move from y = ^/x to y = —y/x. 

If So is a root of A{s, z), then the others will be given by —sq, I/sq, —I/sq. Therefore, we 
can distinguish three cases: 

1. So = —So ^ So = 0, oo ^ 2;^, = oo. 

2. So = 1/ So si = 1 Zc = 1. 

3. So = -I/sq ^ sl = -1 ^ Zc = -1. 

Hence, we have found six Sj, {±1, ±i, 0, oo}, which correspond to three singular points, 
namely Zc = {±1, oo}. By (p2D, we find 

s{z) = ±\/z±V^^^. (34) 

Hence, when z loops around Zc = I 

z-1^ e'^^'Uz-l) s(z) = ±\lz±^/z'^ -l^ ±\l zT^z"^ -1 = (35) 

s[z) 

which is equivalent io K K* . On the other hand, we see that looping around oo may 
correspond to either K —>■ —K or K —>■ K + i7T/2, since 

z e'^'^'z s{z) ~ ±{2z)^^ Ti^z)^^ ~ -s{z). (36) 

Indeed, as z — > oo 

c{z) = cosh2K{z) = ±\/z + l ± Vz^ - 1, (37) 
may have the following asymptotic behaviours 



or 



In the first case we see that 



c{z) ~ ±V2z, 
c{z) ~ ±^^. 

z e^'^'-z =^ c{z) — > -c{z), 



(38) 
(39) 

(40) 



which is equivalent to K ^ K + y, whereas in the other case 

z e^^'z ^ c{z) c{z) (41) 

which is equivalent to K ^ —K. 
Let us set 

H^) = Fugh{K{z)), (42) 

and 

^d(^) = Fiou,{K{z)). (43) 
By ([I^ and (|^), we see that looping around Zc = 1 

T{z) ^ T{z) = Fugh{K{z)) = FMgh{K\z)) = 

Fio^{K{z)) - lnsinh2ir(;2) = J^d{z) - \yis{z). (44) 
In a more symmetric fashion 

J^{z) - i In s{z) ^ ^(5) - ^ In s{z) = - ^ In s{z), (45) 



where we used the fact that s(z) = s{z) ^. Furthermore, by (|T0D and (0), with s 



z 



z ^ e^-z ^ Tn{z) ^ J'niz) + ^. (46) 



4 The Uniformization Equation 

We said above that this geometrical structure is already encoded in Onsager's solution. 
Indeed, the derivative of the free energy, the internal energy per spin U{K) can be expressed 
in terms of elliptic integrals (see for example |TOi) 



UiK) = -^[PfiK)] = -Jcoth2K 



l + k[-}C{h) 

TX 



(47) 



where /C(/ci) is the complete elliptic integral of the first kind 



and 



ki = 4h ^1 = 4-1 kl + {k[f = l. (49) 
+ 1 + 1 
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In terms of F{K{z)) = T{z), using and K = [3J 

U{K)- 



which imphes, after some algebra, that 



a 



+ 



-}C{k 



1), 



(50) 



(51) 



where a = ±1 ^ = z + a\Jz'^ — 1. 



Furthermore, note that \J z — \ IC{ki{z)) is a solution of the uniformization equation for 
the three-punctured sphere A^s = P^(C)\{2; = ±1, oo} 



3 + z' 



4(^2 - ly 



ij{z) 



0. 



(52) 



Therefore, by (|5l|) , J-'{z) solves the following third-order equation 



2z + l 



.2 _ I 



2(^+1)2(^-1) 

Let us denote by H the Poincare upper half plane 



d.J'iz) 



a z + l 
8{z^ - 1)1 



(53) 



H = {TeC\ Imr > 0}. 

The uniformization theorem tells us that the universal covering of the three-punctured sphere 
A^3 is H. Indeed, there exists a holomorphic covering map 



T f-^ Jnir) = z{r), 

such that 

zM = z{t) r(2), 
where r(2) is a subgroup of PSL{2, Z) defined by 

a b 



(54) 
(55) 



r(2) 



c d 



ad — be = 1, a,d = 1 mod 2, 6, c = mod 2 



and 



7r 



ar + b 
CT + d 



(56) 



Hence, = H/T{2). 
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The ratio of two linearly independent solutions of (^) gives the inverse of the uniformiza- 
tion map up to an overall PSL{2, C) transformation 



= Jm\^) = (57) 

The function t{z) is polymorphic: when z loops around one of the punctures, {±l,cx3}, 
iPd{z) and ipiz) undergo a monodromy 



-ipniz) \ _ I 0. b \ f ipoiz) 



(58) 



and correspondingly 

t(z) = , 

ct{z) +d \ c d 

in accordance with (p5|). 



rCz) = ^441^. h r(2). (59) 



5 Potts model 

The Potts model generalizes the Ising model in the sense that each spin Sj can have q 
values, (1, 2, . . . , g), and the Hamiltonian reads 

n = -JY.Ss,sr (60) 

The partition function is given by 

ZiK) = e-^^ = E [K E '^^.^.1 ' (61) 

W} W} \ {id) I 

where as usual K = (33 . Although the problem of finding the exact free energy for the square 
lattice D = 2 Potts model has not been solved, some exact results are known. In particular, 
there exists the analogue of the Kramers- Wannier duality relation , which determines the 
critical temperature 

K, = ln{l + ^). (62) 

For g = 3,4 the phase transition is second order, while for g > 4 the transition is first 
order and the latent heat is known exactly Furthermore, for g > 4 the spontaneous 
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magnetization |T^ and the correlation length at |^0[ are also known. Finally, for g = 3, 4 



the critical exponents are given exactly For a comprehensive review see ||T 
Let us define 

F(K) . Ji,n !!ifl). (63) 

The high temperature expansion reads 

ZuAK) = q'^CiKY'' Y: T{K)\ (64) 

loops 



where 

' = 1' 

q + q — 1 

and / denotes the length of the loop. Conversely, for the low temperature expansion we have 



C{K) = ^'^ \ T{K) = ^, (65) 



ZUK) = qe'^^ Y: e-^'- (66) 

loops 

Note that, in contrast to the Ising model, the length / can be odd. For instance, consider 
g = 3, say blue green and red spins. In the low temperature expansion there are excitations 
corresponding to one green and one red spin next to each other in a background of blue 
spins. In this case / is seven. Therefore, the analogue of (|ll]) is 



K ^ K + 2Txi. (67) 

Under this shift, by (|66D 



Fioy,{K + 27^^) = Fioy,{K) + Ani. (68) 
The generalization of the Kramers-Wannier relation for the Potts model is 

Zhigh{K*) _ Ziow{K) 



2KN 



(69) 



implying that 



where 



FhighiK*) - \nq-2\nC{K*) = Fi^K) - 2K, (70) 



e-^* = TiK) = ' \ - (e^* - l)(e^ - 1) = g. (71) 



Therefore, we see that the critical temperature is indeed given by dl). By virtue of (|7 
we can rewrite Eq. (|70|) in a more symmetrical way, namely 

FMghiK*) - ln(e^* - 1) = Fi^K) - ln(e^ - 1). (72) 
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As before, the mapping K*{K) is involutive. 

The simphcity of transformation (^) makes it easier to find a variable invariant under 
this mapping and (^). We can basically choose either 

w{K) = T{K)T{K*) = e-^e-^" = ^ /"^ ~ ^ (73) 

[e^ + g — 1) 

or 

v{K) =T{K)+T{K*), (74) 
which are related by (|7lD, in particular 

v{K) = {1 - q)w{K) + 1. (75) 

Solving Eg ■ (|73|) we find 

K ^ (1 - g)^ + 1 ± (1 - g)v/(^ - (TT7fF)(^ - (T^) 
^ 2w ■ { ) 

By construction, there are two values of w where the solutions of ( [76| ) coincide. These critical 
values Wc = w{Kc) are given by 

where wi corresponds to the physical singularity and W2 to the complex temperature sin- 
gularity. These are the fixed points of the duality transformation ([7l|) . However, note that 
F{K) is singular as K 0, oo 

K^oo^ Fio^{K) ~ 2K. (78) 

Hence, we argue that there is a further singularity at w{K = 0) = w{K = oo) = 0. 

In order to show the close relationship with the Ising model, we will perform a Mobius 
transformation on w so that the critical points coincide with {±1, oo}. 

The linear fractional transformation that maps wi to 1, W2 to — 1 and w = to oo is 
given by 

Z = <^=^ W = ; r. (79) 

2^w 2^z + {l + q) ^ ^ 

Hence 

[e^-l)2 + g 



z 



(80) 



2^{e^ - 1) ' 
and 

e"" = 1 + ^iz±^/^^). (81) 
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By (|8TD, we see that looping around Zc = 1, K goes to K* 

K{z) = ln(l + ,/qizT Vz^)) ^ (e^(^') - l)(e^(^) _ i) = g ^ K{z) = K*{z), (82) 
where z — 1 = e'^^^{z — 1). Furthermore, choosing the plus sign in (^Tj), as z loops around oo 

z e^'^'z ^ K ^ \n z ^ \nz + 2^1 = K + 27rz, (83) 
that is we retrieve the symmetry 



6 Fisher zeroes 



In []16[, Fisher emphasized the role of the zeroes of the partition function in the complex 
temperature plane in the study of phase transitions. In particular, using Kaufman's expres- 
sion of Z{K) for the two-dimensional Ising model he showed that, in the thermodynamic 
limit, the zeroes are dense on two circles in the tanhi^ or equivalently plane given by 



tanh K =1 + \p2e'\ tanh K = -\ ^ \p2t 



,ie. 



(84) 



These are the antiferromagnetic and ferromagnetic circles respectively and they are related 
by the map K — > —K. 

In terms of z, both circles correspond to the segment joining Zc = —1 to Zc = +1. By 
we have 



2 tanh 



±1 + V2e 



ie 



tanh^ K{e) J 



Then 



which implies that 



AO) 



1 



1 



1 



Therefore, the two circles are constrained on the real axis and 

2 



z{9) 



cos 6 ^ i sin 9 
cos 9 ±i sin 9 , 



+ c.c. 



{\f2 ± COS0 =F % sin( 
((y2±cos^)2 + sin' 



+ c.c. 



(V2±cos0)^-6(v/2± cos6')2 sin 6* + sin 9 
((V2±cos^)2 + sin2^)2 ' 
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implying that 1-2(6') | < 1, as it can be easily checked. Hence, since both 
are in the image, z{^&) maps the interval [0, 27r] into the interval [—1, +1]. 



-1 and Zc = \ 



In |22][23], Martin, Maillard and Rammal conjectured that the Fisher zeroes in the two- 
dimensional g-state Potts model lie on a circle in the complex plane given by 



1 



g — 1 g — 1 

Actually, they conjectured that the locus was given by the condition 



-K* 



-K 



(85) 



(86) 



where K* is the dual of K and K denotes the complex conjugate of K. This yields precisely 
(^5l). Note also that setting q = 2 recovers the ferromagnetic circle of the Ising model. 
Numerical investigations for small lattices at g = 3 and g = 4 |^ provided evidence for 



(pSl). Further progress was made in ||2J], for 3 < g < 8. In p5|, on the basis of numerical 
results on small lattices for g < 10, it was conjectured that for finite lattices with self-dual 
boundary conditions, and for other boundary conditions in the thermodynamic limit, the 
zeroes in the ferromagnetic regime are on the above circle. This conjecture was actually 



proved for infinite g in |26 



In ||2^ , using a general result concerning the partition function zeroes of models displaying 
first order phase transition obtained by Lee |]2^, Kenna proved that the locus for g > 4 is 
indeed given by (p5|) . His argument relies on the fact that the thermodynamic limit and the 



application of the duality transformation ([Til) to the Fisher zeroes should commute. 

We will now show that, in terms of the appropriate modulus z, the locus of the Fisher 
zeroes for the g-state Potts model ( |55[ ) still corresponds to the segment joining Zc = —1 to 

Zc = +1. 

First, note that the locus is equivalent to 



1 + v^e^'. 



(87) 



Indeed, (e 



K\-l 



is a Mobius transformation and thus it maps circles into circles. 



Furthermore, Eq.(p5|) implies that 



yge*^ - 1 



e^(e) _ 1 



(88) 



Finally, by (|8ll) , we immediately recognize that (87) corresponds to z{9) = cos 6 
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z ± Vi^^ = e'^ ^z^- 2ze'^ + e^'^ = z^ z{e) = cos 9. (89) 

Therefore, as in the case of the Ising model, the locus of Fisher zeroes for the Potts model 
corresponds to the segment joining the complex temperature singularity Zc = —1 to the 
physical singularity Zc = +1. 



7 Conclusions 

In summary, we have seen that the introduction of a proper modulus or uniformizing coor 
dinate z, which is given by 

1 + sinh^(2ir) 



in the Ising case, and 



2sinh2(2ir) 
(e^ - 1)2 + g 



(90) 



(91) 



2^{e^ - 1) ' 

in the Potts case, unveils an intriguing correspondence between the D = 2 Ising and Potts 
models. Indeed, their moduli spaces are both equivalent to the three-punctured sphere 

= pi(C)\{^ = ±l,oo}. (92) 

The moduli (^) and (^) are characterized by the invariance under the symmetry and 
duality transformations (0) (|llD (0) and (|67D ([7TD respectively. 



This connection is further strengthened by the fact that the locus of Fisher zeroes in 
both models is given by the segment joining Zc = —1 to Zc = +1. 
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